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This paper investigates the use of the inverse first-order reliability method (inverse- 
FORM) to quantify the uncertainty in the remaining useful life (RUL) of aerospace compo- 
nents. The prediction of remaining useful life is an integral part of system health prognosis, 
and directly helps in online health monitoring and decision-making. However, the predic- 
tion of remaining useful life is affected by several sources of uncertainty, and therefore it is 
necessary to quantify the uncertainty in the remaining useful life prediction. While system 
parameter uncertainty and physical variability can be easily included in inverse-FORM, 
this paper extends the methodology to include: (1) future loading uncertainty, (2) process 
noise; and (3) uncertainty in the state estimate. The inverse-FORM method has been used 
in this paper to (1) quickly obtain probability bounds on the remaining useful life predic- 
tion; and (2) calculate the entire probability distribution of remaining useful life prediction, 
and the results are verified against Monte Carlo sampling. The proposed methodology is 
illustrated using a numerical example. 

I. Introduction 

The use of accurate system health management tools is indispensable for safety-critical and mission- 
critical aerospace systems. It is important to constantly monitor the performance of these systems, identify 
faults (diagnosis), predict possible failures in the near future, and quantify the remaining useful life (progno- 
sis) in order to aid online decision-making. Sometimes, it may be challenging to perform health monitoring 
on the whole system due to its sheer complexity, and therefore, diagnosis and prognosis need to performed 
on individual components which constitute the overall system. In this approach, mathematical models are 
developed for individual components, and then the component models are integrated to form the overall 
system. These models can be used in the health monitoring to guide in model-based diagnostic^ and prog- 
nostics . 2 Alternatively, data-driven approached are also available for health monitoring, but they are not 
discussed in the present paper. 

Uncertainty quantification and management are important aspects of health monitoring, especially due to 
the presence of several unknown factors which affect the performance of the system of interest. It is not only 
important to develop robust algorithms for diagnosis and prognosis, i.e. , accurately perform diagnosis and 
prognosis in the presence of uncertainty, but also important to quantify the amount of confidence in the results 
of diagnosis and prognosis to enable in-flight decision-making capabilities. Sankararaman and MahadevaiP® 
developed statistical (both frequentist and Bayesian) approaches to quantify the uncertainty in the three 
steps of diagnosis (detection, isolation, quantification/estimation) in an online health monitoring framework. 
There have also been a few papers® 7 8 9 which discuss uncertainty propagation in prognosis; however, many 
of these papers are either suitable only for offline prognosis or they do not provide a comprehensive treatment 
of uncertainty .® 8 Further research is necessary to develop rigorous uncertainty quantification methods which 
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are suitable for online prognostics and decision-making. The rest of the paper focuses on one particular aspect 
of prognostics, i.e., quantifying the uncertainty in the remaining useful life (RUL) prediction. Since state 
space equations are suitable for representing time-dependent behavior, the governing differential equations 
of the component (of the system of interest) are transformed into equivalent state space representations that 
collectively serve as the component model for performance prediction and RUL computation. The first and 
foremost challenge with respect to the estimation of RUL uncertainty is the characterization of uncertainty 
in the future loading. Though a few publications have addressed in this issuefSE 1 further research is 
needed in this direction. However, the focus of the present paper is not on uncertainty characterization, 
and therefore it is assumed that future loading variability has been characterized based on already existing 
loading profiles and data sets. The other sources of uncertainty include the process noise of the state space 
model (a component of modeling error), the uncertainty in the current state estimates, uncertainty in the 
model parameters, etc. Note that state estimation is a prerequisite to prognosis, and a filtering algorithm 
is used along with measurement data to estimate the current state. The data uncertainty associated with 
these measurements directly contributes to the estimated uncertainty in the component states. 

Conventionally, prognosis has been performed using sampling-based approaches (for example, particle 
filtering without the correction-step); since such methods may require several samples or particles in order 
to predict the RUL, and may not be suitable for online decision-making. Moreover, sampling-based based 
algorithms are “non-deterministic” ; for example, every instance of a particle filter may lead to a different 
probability distribution for the RUL. On the other analytical methods such as the first-order reliability 
method (FORM), though approximate, can overcome both of these challenges. In particular, the inverse 
first-order reliability method (inverse-FORM) can be used to calculate probability bounds on RUL (in 
fact, the entire probability distribution of RUL at discrete probability levels), which are useful for decision- 
making. Though such analytical methods have been used for time-dependent reliability calculation , 12 13 
their application to state space models poses several new challenges. The primary contribution of this 
paper is to extend the inverse FORM approach to include state space models and compute the uncertainty 
in the RUL prediction. The inverse-FORM approach not only requires fewer prognostic evaluations in 
comparison with the sampling-based approaches but can also produce repeatable calculations, i.e., the exact 
same probability distribution on every repetition of the algorithm. While the former directly aids in online 
prognosis since fewer evaluations would lead to quicker calculations, it is worth noting that the latter feature 
is an important criterion for existing verification, validation, and certification protocols in the aerospace 
domain. Therefore, investigating such analytical algorithms allows us to move a step closer towards adopting 
prediction algorithms, by meeting the needs of the current certification process. The rest of this paper is 
organized as follows. Section [TT] discusses the various challenges involved with using inverse FORM along 
with state space models, and develops computational methods to address them so that the uncertainty in the 
RUL prediction can be computed. Section IIIII illustrates the proposed methods using a numerical example 
and conclusions are presented in Section IIVI 

II. Quantifying Uncertainty in RUL Prediction 

First, the inverse-FORM approach is briefly explained in subsection IHIAI In subsection llllBl a generalized 
framework for the computation of RUL is described, and then, the challenges of using FORM with such state- 
space models is explained in subsection IHICI The proposed method is described in subsection IHIDI and finally, 
the inverse-FORM algorithm is summarized in subsection IIIIEI 

A. The Inverse FORM Approach 

Consider a generic computational model Y = H(X), which is used to represent the performance of an 
engineering system. The input is a vector and hence denoted in bold as X , whereas the output Y is a scalar. 
The model G is deterministic, i.e. for a given realization of X , there is a corresponding output, which is 
a realization of Y . The inputs X are uncertain, and this leads to uncertainty in the output Y . A generic 
realization of X is denoted as x , and a generic realization of Y is denoted as y. Given the probability 
distributions of X, FORM and inverse-FORM are tools for uncertainty propagation, i.e., to calculate the 
cumulative distribution (CDF, denoted by Fy(y)) or the probability density function (PDF, denoted by 

The focus of inverse-FORM'Sl is to select the value of y c such that the cumulative distribution function 
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Figure 1: Estimating MPP in inverse-FORM 


(CDF) value measured at y c is equal to a given probability value a, i.e. select y c such that Fy(Y = y c ) = a. 
The basic concept is to “linearize” the model H so that the the output Y can be expressed as a linear 
combination of the random variables. First, all the random variables are transformed into uncorrelated 
standard normal space; in this transformed space, the point of linearization is chosen such that it is closest to 
the origin, implying that it has the highest probability of occurrence (Most Probable Point or MPP). Further, 
in the transformed space, the output Y is also a normal variable because a linear combination of normal 
variables is normal. Therefore, the CDF value Fy(Y = y c ) can be computed using the standard normal 
distribution function in terms of the minimum distance (denoted by /3 in Fig. [Tj) as Fy(Y = y c ) = $(— /3). 
The transformation of random variables X into uncorrelated standard normal space ( U ) is denoted by 
U = T(X) 1 and the details of the transformation can be found in Haidar and MahadevanfU An iterative 
procedure (in each iteration, (1) the transformation is repeated; and (2) the point of “linearization” changes 
when the U H” is non-linear and X are non-normal variables) for inverse-FORM can be found several 
publications (for e.g., Li and FoschP^). The iterative procedure starts at an arbitrary x and continues until 
all the optimality conditions are satisfied. 

While inverse-FORM has commonly been used to calculate y c corresponding to a given a, it can also be 
used to compute the entire CDF by repeating the inverse-FORM analysis for multiple values of a (different 
probability levels). For example, if inverse-FORM is performed at 10 different values of a, the corresponding 
y c values are calculated, and an interpolation scheme can be used to calculate the entire CDF (which can 
be differentiated to obtain the PDF). 

B. Generalized State Space Model for RUL Prediction 

State space models are generally suitable to represent time-dependent behavior, and therefore, such models 
are commonly used for prognosis and RUL prediction. A general state space model can be represented as: 

x(t) = f (t, x(t), G(t), u(i), v(i)) (1) 

where x(t) £ M™ x is the state vector, Q(t ) £ R™ 8 is the parameter vector, u(t) £ K™“ is the input vector, 
v(t) £ R n ’' is the process noise vector, f is the state equation. 

The goal of prognostics is to predict the remaining useful life (RUL) at any generic prediction time 
(tp). The first step is to estimate the state at this prediction time, i.e., x(tp); such estimation is typically 
performed using filtering approaches (particle filtering, Kalman filtering) using the data available until the 
prediction time (tp). Let y (t) £ R n “, n(i) £ I” n , and h denote the output vector, measurement noise 
vector, and output equation respectively. Then,, 

y(t) = h (t, x.(t),0(t), u (t), n (t)) (2) 

Note that the output equation or data is not used in the prognosis stage, since the focus is on predicting the 
future and the associated uncertainty. 
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Prognostics and RUL prediction is concerned with the performance of the component that lies outside a 
given region of acceptable behavior. The desired performance is expressed through a set of n c constraints, 
Ceol = where Ci : R™* x R™ 8 x R™“ — >• B maps a given point in the joint state-parameter space 

given the current inputs, (x(t), 9(t),u(t)), to the Boolean domain B = [0, 1], where a(x.(t),9(t), u(t)) = 1 if 
the constraint is satisfied, and 0 otherwise. 

These individual constraints may be combined into a single threshold function Teol '■ R"* x R” 8 x R” u — >• 
B, defined as: 


T EO L{x-{t),0{t),u(t)) 


i, oe{cj(x(t),0(t),u(f))gi 

0, otherwise. 


( 3 ) 


Teol is equal to 1 when any of the constraints are violated. Then, the End of Life (EOL, denoted by E ) 
at any time instant tp is then defined as the earliest time point at which this occurs: 


E(t P ) = inf {t e R : t > t P A T EOL {x(t), 0(t),u(t)) = 1} 


( 4 ) 


The Remaining Useful Life (RUL, denoted by R) at time instant t P is expressed as: 

R(tp ) = E(tp) — t p . 


( 5 ) 


Thus, it is clear that R(t P ) depends on (1) present time (tp); (2) present state estimate (x(t P )); (3) future 
loading (u (t P : t P + R(t P ))); (4) parameters (B(t P : t P + R(t P ))); and (5) process noise (v(t P : t P + R(t P ))). 
Since these quantities may be uncertain, RUL also becomes uncertain, and as described earlier, the goal 
of this paper is to quantify this uncertainty. For the sake of illustration, assume that the parameters are 
known deterministically; even if they are not, parameter uncertainty can be treated to the other sources of 
uncertainty. 

Note: The symbol “u(tp : tp + R(tp ))” refers to the entire loading from time “t = tp” until time 
“t = tp + R(t P )” . Similarly, the parameter values and process noise values during this time period can be 
written as “6(tp : t P + R(tp))” and “v(tp : tp + R(tp))” respectively.. 


C. Challenges 

The first challenge is that the dependence of R(tp) on other uncertain quantities is implicit, i.e., the RUL 
depends on the process noise and loading, and in turn, “how long the process noise and the future loading 
need to be considered during prognosis” directly depends on the RUL itself. 

In the past time-variant random loading has often be decomposed into two parts: (1) time-variant 
deterministic component; and (2) time-invariant random component. Therefore, u (tp : tp + R(t P ))) can 
be replaced with u E (the equivalent, time-invariant random components), and thereby, such decomposition 
facilitates the treatment of future loading uncertainty. However, the treatment of process noise (usually, but 
not necessarily, considered to be normally distributed, i.e., v(f) ~ N( 0, E^)) is not trivial, because v(U) and 
v(< 2 ) need to be treated as two different random quantities (there may be correlation across time, which is 
not considered in this paper). So, there are “R(tp)x time-discretization” number of process noise terms; 
note the presence of implicitness since this is again a function of the RUL itself. 

The next challenge is to identify the function “H” discussed earlier in Section [A] The obvious choice 
is to use the state evolution equation, i.e, Eq. |T] as “H” . Then, this would mean performing uncertainty 
quantification across every time-step, thereby predicting the entire probability distribution of all the state 
variables (x(t P : tp + R(t P ))) continuously as a function of time. However, this is a cumbersome process 
because of the fact that there may be several states which are highly statistically dependent (sometimes, they 
may be even functionally dependent), and therefore, it is important to retain the correlation information all 
along. Further, these state variables may not be normally distributed, and every inverse-FORM analysis 
would require a transform from correlated non-normal space to uncorrelated standard normal space, thus 
complicating the procedure. In the next subsection, a statistical method is proposed in order to overcome 
these challenges and facilitate the implementation of inverse-FORM for RUL calculation. 


D. Proposed Method 

Since the focus of this paper is to compute the uncertainty in the RUL, it is proposed that a new deterministic 
transfer function “G” be constructed, so that the inputs are realizations of uncertain quantities, and the 
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output is the corresponding realization of RUL. This function can be written as: 

R(t P ) = G(x(t P ), u B , v(f P : t P + R(t P ))) (6) 

This function can be evaluated by continuously solving Eq. [T] until failure occurs (when Teol changes the 
value from 0 to 1). While this “G” is similar to “iJ” in Sec [A] there are two more challenges that need 
to be addressed. First, the above equation is still implicit, because the number of process noise terms is 
dependent on R(tp). Second, if R(tp) is discretized into several hundreds of steps, then there are several 
hundreds of inputs to “G” . This poses a computational challenge to inverse-FORM because of the need to 
compute the gradient that may now have several hundreds of terms. Therefore, to obtain the entire PDF 
(repeating inverse-FORM for 10 different probability levels), several thousands of evaluations of “G” may 
be necessary. Therefore, this approach becomes more expensive than sampling methods, and not suitable 
for online health monitoring. 

This paper proposes a new likelihood-based method to overcome this challenge. Suppose that the same 
realization of process noise v(t) = v E is considered for all time-steps t > tp, i.e. , the process noise is time- 
invariant (constant at all time steps after prediction time tp). Then, the likelihood for the occurrence of this 
event can be written as: 

t=tp+R(tp) 

L(v E \R(t P )) oc /v(t)(V(t) = v E ) (7) 

t=tp 

where /v(t)( v (0) is the probability density function of V(f). Also note that the likelihood function is 
conditioned on the RUL and written as R(tp). Further, the above equation assumes that the process noise 
values at two different times are independent of each other. If any statistical dependence is unknown, then it 
can be easily included in the above equation by conditioning appropriately. Having calculate the likelihood 
of this event, the PDF of v E can be calculated asP^ 


fvE{y E \R(t P )) 


L(v E \R(tp)) 
f D L(v E \R(tp))dv E 


( 8 ) 


In Eq. [HI the domain of integration D is chosen such that v E € D if and only if L(v E ) ^ 0. 

Now, suppose that, in Eq. [Gl the time- variant process is replaced with a time-invariant constant. In other 
words, v(tp : tp + i?(£p))) is replaced with v E , as: 

R{tp) = G(x(t P ), u £ , v B ) (9) 

The above discussion means that a suitable probability distribution needs to be chosen for v B so that the 
effect of propagating this distribution through Eq. [9] is equivalent to propagating the original distribution of 
V(t) through Eq|G] In fact, Eq.[7]and[H]accomplish this goal by matching corresponding probability densities. 
Hereon, v E is referred to as equivalent time-invariant process noise and the distribution corresponding to 
the PDF f V E{y E \R(tp )) is referred to as the equivalent invariant process noise distribution. 

Now the function “G” is not implicit, because v E is simply a time- invariant constant and R(tp) does not 
appear on the right hand side of Eq. [9] On the other hand, the PDF of v E is conditionally dependent on 
R(tp), and this dependence can be expressed explicitly. In other words, for a given realization of v B , Eq. [9] 
can be used to compute R(tp), which can then be used to calculate the PDF f- V E(v E \R(tp)). Hence, any 
iterative algorithm for inverse-FORM can now be applied for RUL estimation by performing one additional 
step, i.e., calculating the probability density function of the current value of v B , which is then used to 
transforming the variables to the standard normal space. 


E. Algorithm: Summary 

The section summarizes the algorithm used in this paper and explains how inverse-FORM is to compute the 
value of RUL that corresponds to a given probability level, i.e., compute r for a given a such that Fp(r) = a. 
To begin with, calculate j3 such that a = Phi(—/3), and use the following algorithm to calculate r. 

1. Start with an initial guess of all the quantities which RUL depends on, as in Eq. [9] Initialize j = 0, 
and let the vector of initial guesses be denoted by xP 
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2. Transform x J into standard normal space w- 7 , using any of the commonly used transforms.^ 5 The 
probability distribution of v E needs to be calculated as in Eq. [8l Note that the value of R(tp) is also 
calculated for each iteration and this facilitates the computation in Eq. [8] 

Calculate the gradient (a) in the standard normal space, and update the location using the gradient, 
i.e., « 1+1 = jfj x p. 

Re-transform into original space using the inverse of the previously used standard normal transforma- 
tion, and obtain x^ +1 . 

5. Repeat steps 2 — 4 until convergence, i.e., x = x^ +1 and \x^\ = (3. 

The following section presents a numerical example to illustrate the above algorithm; first the state-space 
model is discussed in detail, and then the proposed method is used to quantify the uncertainty in remaining 
useful life prediction. 


3. 

4. 


III. Numerical Example 

The numerical example consists of a lithium-ion battery, which is an important component of the rover 
test-beef 1 ® being studied at NASA Ames Research Center. Though the method is illustrated using a battery 
model, it is general enough to be applicable to state space models in several engineering domains. 

A. Description of the Model 

The battery model is based on an electrical circuit equivalent shown in Fig. [2] similar to models presented by 
Chen and Rincon- MoraP 0 The large capacitance Cb holds the charge qb of the battery. The Rcp-Ccp pair 
captures the major nonlinear voltage drop due to concentration polarization, R s captures the so-called I-R 
drop, and R p models the parasitic resistance that accounts for self-discharge. This simple battery modePSIH 
is enough to capture the major dynamics of the battery, but ignores temperature effects and other minor 
battery processes. 


R CP 



The state-of-charge, SOC , is computed as 


soc = 1 - qmax qb ( 10 ) 

max 

where qb is the current charge in the battery (related to Cb), q m ax is the maximum possible charge, and 
Cmax is the maximum possible battery capacity (i.e., nominally, its rated capacity). The concentration 
polarization resistance is a nonlinear function of SOC: 

Rcp = Rcpo + Rcpi exp(J?cP2(l - SOC)) (11) 

where Rcpo, Rcpi, and Rcp i are empirical parameters. The resistance, and, hence, the voltage drop, 
increases exponentially as SOC decreases.®! 
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Voltage drops across the individual circuit elements are given by 


Vb = q b /C b (12) 

Vcp = qcp/Ccp (13) 

V p = V b - V CP (14) 

where qcp is the charge associated with the capacitance Cc p- The terminal voltage of the battery is 

V = V b ~V CP - R s i (15) 


where i is the battery current at the terminals. Currents associated with the individual circuit elements are 
given by 


ip = V p /R p (16) 

i b = i p + i (17) 

icp = ib — Vcp/Rcp (18) 

The charges are then governed by 

q b = ~i b (19) 

qcp = icp ( 20 ) 


It is of interest to predict the end-of-discharge as defined by a voltage threshold Veod • So, Ceol consists 
of only one constraint: 

ci : V > Veod (21) 

The parameters of the battery model are assumed to be deterministic and are shown in Table [Tj 


Table 1: Battery Model Parameters 


Parameter 

Value 

c b 

9844 

Rs 

0.143014 

R P 

500 

Ccp 

70.3767 

Rcpo 

0.019829 

Rcp i 

3.68606 x 10" 14 

RCP2 

31.9213 

Qmax 

41400 

Cmax 

6900 


B. Verification of the Equivalent Time-invariant Approach 

Before quantifying the uncertainty in remaining useful life (RUL) prediction, this section focuses on verifying 
the proposed equivalent time-invariant approach, by suppressing all the other sources of uncertainty except 
process noise. Since there are two state variables in this example, two process noise terms are present. The 
two initial states are assumed to [41. 4x 10 -3 , 0], and time- invariant PDFs are chosen for the two process noise 
terms with normal distributions of N( 0, 1) and 1V(0,0.1) respectively; these are the probability distributions 
at any time instant and used to calculate the probability distributions in Eq[8] Even though small values of 
standard deviation are chosen, they have a significant impact on the RUL distribution, as shown below. 

The statistical correctness of the proposed equivalent time-invariant approach is independent of the 
uncertainty propagation technique (sampling-based or analytical) used. Therefore, for the sake of rigorous 
verification, Monte Carlo sampling (MCS) is used to compare the CDFs obtained by using the original 
process noise distribution and the equivalent invariant process noise distribution, i.e., 

1. Proposed time-invariant approach, i.e. same value of process noise at all time steps, as in Eq. [9] The 
resultant CDF is indicated as “Invariant Noise” in Fig. [3j 
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Remaining Useful Life 

Figure 3: Verification of the Proposed Equivalent Time Invariant Approach 


2. Actual process noise, i.e., sampling independent realizations of process noise in each and every time 
step with Eq. [G] In this case, 3 -a bounds of the Monte carlo estimate (due to the use of finite number 
of samples) is also computed. The resultant CDFs are indicated as “Actual” in Fig. [3] 

Note that both the above methods are sampling-based; while the former uses the same process noise values 
at all time steps, the latter selects a random sample of the process noise in each time step. The comparison 
between the above two is shown in Fig. [3j and it is seen that the proposed equivalent time-invariant approach 
compares well with that obtained using the true process noise. Thus, the use of time-invariant process noise 
has been validated, and this can now be used in inverse-FORM for quantifying the uncertainty in RUL. 

C. Quantification of Uncertainty in RUL 

Now that the proposed likelihood-based method has been verified, the next task is to replace Monte Carlo 
sampling with inverse-FORM. This is now easy because the original inverse-FORM can be used, since all 
the implicit relations have been eliminated. The different sources of uncertainty considered in this example 
are described below: 

1. State Uncertainty: Typically, the state estimation, which is an inverse problem, is addressed using 

a filtering technique which can continuously estimate the uncertainty in the state when measurements 
are continuously available as a function of time. In this paper, the state estimation is not explicitly 
carried out. The state values are assumed to be available, and the uncertainty in the states are 
predetermined based on the authors’ past experiences with the use of filtering techniques for the above 
described problem. There are two state variables (charge in the battery and charge associated with the 
capacitance) in this example and at any time instant, they are assumed to be normally distributed with 
a specified mean; for example, the mean of the initial states are set as [4.14 x 10 4 , 0]. For the purpose 
of illustration, three different values of CoV (Coefficient of variation, defined as the ratio between 
standard deviation and mean) are considered - 0.01, 0.02, 0.03 — and the analysis is repeated for 

each CoV value. 

2. Loading Uncertainty: In this example, a constant amplitude loading is considered for the pur- 
pose of illustration. The amplitude is considered to be normally distributed (IV(1.375, 1/6)), and this 
distribution is truncated at a specified lower bound (0.75) and an upper bound (2.00). 

3. Process Noise: In each time step, the two process noise values are assumed to normally distributed 
as IV(0, 1) and IV(0,0.1) respectively. 

The proposed inverse-FORM approach (including time-invariant process noise) is used to quantify the 
uncertainty in RUL and calculate the entire CDF of RUL (for each CoV corresponding to state uncertainty) 
continuously as a function of time. The 98% probability bounds of RUL for three cases (CoV of state 
uncertainty = 1%, 2%, 3%) are shown in Figs. ITal ITcl 
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(b) 98% Bounds for CoV=2% 




Figure 4: Probability Distribution of RUL 


The results from inverse-FORM also seem to be comparable against the Monte Carlo solution; for exam- 
ple, the CDFs from the proposed approach and the Monte Carlo simulation approach for RUL calculation 
at initial time and 3% CoV are shown in Fig. I4dl In fact, a similar agreement was confirmed at multiple 
time instants for all the three CoV values, thereby verifying the proposed approach. 

IV. Conclusion 

This paper discussed the use of the inverse first-order reliability method (inverse FORM) in order to 
compute the remaining useful life of aerospace components. The computation of remaining useful life is 
important in the context of component-level and system-level prognosis, and hence directly useful for online 
operations and decision-making. Conventionally, sampling-based algorithms have been used for quantifying 
the uncertainty in prognostic calculations, and may require several thousands of system evaluations in order 
to quantify the uncertainty in remaining useful life (RUL) with reasonable accuracy. On the contrary, 
the inverse FORM approach is an analytical algorithm and computes the uncertainty using a few system 
evaluations, and therefore is suitable for online prognosis. Further, the inverse FORM approach is invariant 
on repetition (as against sampling methods like Monte Carlo analysis), and hence may be preferred for 
system validation and certification purposes. 

While this paper addressed a few challenges involved in using inverse-FORM for prognostics and calcu- 
lating the uncertainty in RUL prediction, there are several other issues that need to be addressed. Practical 
systems are commonly subjected to different types of variable amplitude loading profiles such as block load- 
ing, Markov processes, general random processes, etc., and therefore, the proposed methods for uncertainty 
quantification need to be extended to consider variable amplitude loading. The assumption of constant am- 
plitude loading implies that loading uncertainty is described using a single random variable, whereas variable 
amplitude loading profiles need to be described using multiple random variables, which not only increases 
the dimensionality of the problem, but also affects the uncertainty in the RUL prediction. Therefore, the 
impact of including variable amplitude loading on the uncertainty bounds of the RUL prediction needs to 
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be investigated. Further, sensitivity analysis needs to be performed so that the contributions of the different 
sources of uncertainty to the overall uncertainty in RUL can be quantified. This paper did not consider 
the effect of model form uncertainty on prognosis; future research needs to quantify model form uncertainty 
and develop a method to rigorously account for model uncertainty in prognosis and RUL calculations. Fi- 
nally, it is also necessary to quantify the robustness of the proposed approach, and thereby investigate the 
applicability of the methodology to practical engineering systems. 
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